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(54) Precede de f iltrage spatial du bruit dans une image numerique, et dispositif mettant en 
oeuvre ce precede 


(57) L'invention concerne un precede de filtrage 
spatial du bruit dans une image numeriqua Ce precede 
comprend ia localisation d'un pixel courant [Po(xo. Yo)] 
d'intensite bruitee [lo(xo.yo)l : la definition d'une fenetre 
glissante [WD{xo. yo)] centres sur le pixel courant ; la 
seiection d'un nonnbre (N). fixe ou aleatoire, de filtres 
el6mentaires aleatoires orientes [F^, Fa,,.. F^) ayant un 

axe d'orientation [PqZi. P0Z2 ^o^h] passant par ce 

pixel courant, et munis de caracteristiques dont au 
moins une est aleatoire, incluant : un angle d'orientation 
[©T. 02- ®n] raxe d'orientation et un axe de 
reference (OX), et une etendue de support definie en 
nonnbre de pixels ; la combinaison des intensites des 
pixels sur les supports des fittres elementaires indivi- 
duellement pour fournir des intensites elementaires ; et 
la combinaison des intensites elementaires pour fournir 
une intensite fiitree [Ro(xo. Vo)) a appliquer au pixel 
courant. 

L'invention concerne aussi un dispositif de traite- 
ment d'image mettant en oeuvre ce procede. 
Aooiication : Imagerie medicale numerisee 
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Description 

Unvention conceme un precede de fiitrage spatial 
du bruit dans une image numerique. 

^^nverttion concerue egalefpent un dispositif met- 5 
tant 60 oeuvre ce precede. 

L'invention'trouve particulierement son application 
dans le traitement dimages medicales formees dans un 
orode fluoroscopique de rayons X, au moyen d'un sys- 
t^e cu le faisceau de rayons X a une faibie intenstte. w 
fournissant par suite des images bruitees et peu con- 
trastees qui doivent etre filtrees afin d'en extraire le 
bruit* sans en deteriorer les details. 

Le mode fluoroscopique de rayons X est utilise par 
exemple pour suivre en temps reel une operation medi- is 
cale durarrt laqueile un outil de d[am§tre extrSmement 
petit tel qu'un catheter est introduit ou deplace darDs la 
zone observee sur un patient. 

Un precede de fiitrage du bruit dans une image 
numerique est dej^ connu par la. publication intitulee "A 20 
new class of Detail-Preserving Filters for Image Proces- 
sing" par Ari NIEMINEN. et alii, in "IEEE TRANSAC- 
TIONS ON PATTERN ANALYSIS AND MACHINE 
INTELLIGENCE VOLPAMI-9. N«1. JANUARY 1987, 
PR74-90". 25 

Cette publication d^rit un precede de fiitrage utili- 
sant une nouvelle classe de filtres spatials. Cette nou- 
velle classe comprend des structures de filtres lineaires 
a reponse impulsionnelle, appeles "FiR". Le precede 
divulgue par ce document cite inciut des combinaisons 3o 
des resultats des filtres lineaires dans des filtres 
medians, avec plusieurs niveaux de combinaison. II est 
enseign6 qu'un tei precede permet de preserver les 
bords des objets dans I'image, et ainsi de preserver les 
details. Ce precede inciut la definition d'une fenetre glis- 35 
sarrte sur I'lmage numerique k traiter centree sur un 
pixel courant. et la definition de masques dans cette 
fenetre pour permettre la realisation d'un certain nom- 
bre de moyennes des intensites des pixels situes autour 
du pixel central de la feneti-e. Par exemple. une 40 
moyenne est realisee sur les intensites des pixels situes 
sur la meme ligne de pixels que le pixel central, une 
moyenne est realisee sur les intensites des pixels situes 
sur la m^e colonne que le pixel central ; les moyennes 
d'intensite peuvent aussi Sfre realisees sur des pixels 
des diagonaies de la fenSfre, cette derniere etant can-ee 
dans I'exemple decrit Les resultats de ces moyennes. 
qui constituent des filtrages lineaires (FIR) appeles fil- 
trages de base, sont combines dans des filtres medians 
(FMH), Par exemple. les moyennes calcul^es sur les so 
intensites des pixels de la ligne du pixel central a gau- 
che et a droite du pixel central respectivement sont 
combinees dans un filtr e median avec I'intensite du pixel 
central, et un premier resultat est fourni ; puis les 
moyennes caiculees sur les pixels de la colonne du ss 
pixel central vers !e haut et vers le bas respectivement 
sont combinees avec I'intensite du pixel cerrtral dans un 
fittre median, et un second resultat est fourni ; puis les 
premier et second resultats sont combines dans un fiitre 


median avec I'intensite du pixel central, et un autre 
resultat est fourni. Selon le nombre de filtrages de base 
calojles dans la fen§tre. et le nombre d'entrees des fil- 
tres medians, le nombre d'^tapes de combinaisons pour 
o}:tenir le resultat appeiees niveaux de combinaisons. 
peut etre plus ou moins grand. 

Un auti'e precede de fiitrage du bruit dans une 
image numerique est 6galement dej^ connu par !a 
publication intitulee "Performance Evaluation of Some 
Noise Reduction Methods" par Wen-Yen WU et alii, 
dans "GRAPHICAL MODELS AND IMAGE PROCES- 
SINGS. VOL54. N^a. March 1992. pp. 134- 146. 

Cette seconde publication d6crtt differents prece- 
des de fiitrage du bruit incluant d'une maniere generale 
d'abord la selection, selon des entires donnes. de 
pixels dans une image numerique sur I'intensite des- 
quels on realise une combinaison. Cette combinaison 
peut §tre une combinaison iineaire telle qu'une 
moyenne. A cet effet, ce precede inciut la definition 
d'une feneti-e glissante sur I'image numerique k traiter. 
centree sur un pixel courant, et la selection de pixels 
dans cette fen§tre. Un critere de selection peut porter 
sur I'intensite des pixels selectionnes, Dans un exemple 
(SIG). les pixels selectionnes sent ceux dent la diffe- 
rence d'intensite vis-^-vis du pixel courant situe au cen- 
tre de la fenetre, est inferieure a un seuil donne qui 
depend de la variance du bruit (o^) ; dans un autre 
exemple (KNN). un nombre K donne de pixels ayant les 
intensites les plus proches de celtes du pixel courant est 
selectionne. Apres I'operation de selection des pixels 
dans la fen§tre glissante carree, ieur intensrte est com- 
binee soit en effectuant une moyenne (MEAN), sort par 
un fiitrage MEDIAN (MED). 

Un probl6me qui se pose dans le traitement du bruit 
des images en fluoroscopie est que ces images con- 
tiennent en general des details ti-es fins ou la represen- 
tation d'objets tres petits. comme par exemple un 
catheter introduit lers de I'operation par le praticien, ce 
catheter ayant la forme d'un fil dont la dimension trans- 
versale est de quelques pixels. Dans ce domaine, it est 
done indispensable d'avoir recours k une method e de 
fiitrage du bruit qui conserve tes details. 

Toutes les metiiodes decrites dans les deux publi- 
cations citees precedemment au titre d'etat de la techni- 
que sont capables de fournir un fiitrage de I'image 
bruitee en conservant les details. 

Cependant, les methodes connues des documents 
dtes posent toutes un meme probl^me ; les image's 
reconstruites k I'aide de ces metiiodes monti-ent des 
TACHES (PATCHES). Ceci est du au fait que ces 
methodes reposent sur I'utilisatien d'une fenetre glis- 
sante dont le centre est deplace d'un pixel sur llmage 
apres la mise en oeuvre de chaque ensemble d'etapes 
de calcul decrit dans Tune ou I'autre des publications 
cities. Ainsi. I'aire de recouvrement des fenetres cor- 
respondant a deux pixels adjacents est important e. d'ou 
il resulte que la contribution du bruit residue! dans ces 
deux feneti-es est correlee. 

Un but de la presente invention est de proposer un 
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procede de filtrage spatial qui ne gen^e pas de 
TACHESde bruit. 

Un but particulier de llnvention est de proposer un 
tel procede capable de produire une elimination tres 
pouss^e du bruit* avec preservation des contours des 
objets dans rinr^age, et avec preservation du contraste 
de rimage. 

Ces buts sont atteints par un procede de filtrage 
spatial du bruit dans une image ayant la forme d*une 
matrice bicOmensionnelle de pixels ayant des valeurs 
d'intensite numeris6es bruitees. ce procede compre- 
nant des etapes de : 

a) selection, dans une fendtre glissante centree sur 
un pixel courant d'un nombre (N ^ 0) de f litres ele- 
mental res orientes al^atoires ayant des supports 
avec un axe d'orientation passant par le pixel cou- 
rant et des caracteristiques dont une ou plusieurs 
ont des valeurs aleatoires. 

b) combinaison des intensites des pixels sur lesdits 
supports individuellement. fournissant une intensite 
el6mentaire. et combinaison des intensites elemen- 
taires, fournissant une intensite filtree a appliquer 
au pixel courant, 

Ce procede presente I'avantage de ne pas generer 
de TACHES de bruit residue! dans les images, et de ne 
pas faire apparaTtre non plus des motifs dus a la forme 
des supports des f litres. 

Un dispositrf pour Taide k I'observation d'une image 
medicale comprend : 

un systeme pour fournir des donnees d'une image 
sous forme de matrice bidimensionnelle de pixel 
ayant des valeurs dintensite numerisees, 
un systeme d'affichage pour i'affichage de I'image, 
un systeme de traitement d'images ayant acces 
aux donnees d'image. et au systeme d'affichage. 
ce systeme de traitement d'image etant arrange 
pour effectuer au moins des etapes de : 

a) selection, dans une fen§tre glissante cen- 
tree sur un pixel courant d'un nombre (N ^ 0) 
de filtres elementaires orientes aleatoires 
ayant des supports avec un axe d'orientation 
passant par le pixel courant et des caracteristi- 
ques dont une ou plusieurs ont des valeurs 
aleatoires. 

b) combinaison des intensites des pixels sur 
lesdits supports individuellement fournissant 
une intensite elementaire, et combinaison des 
intensites elementaires, fournissant une inten- 
site filtree k appliquer au pixel courant 

Ce dispositif presente I'avantage de restaur er une 
image medicale bruitee et peu contrastee a I'origine. en 
preservant les plus fins details, et sans TACHES de 
bruit. 

Llnvention est decrite d-apres en detail en refe- 


rence avec les figures schematiques annexees dont : 

la FiG.1 represente un dispositrf a rayons X pour la 

formation d'images radiologiques ou fluorosccpi- 
5 ' ques, 

la F1G.2 represente le passage d'une. JenStre 

WD(xy) dans une image J bruitee. 

les F1G.3A ^ 3D montrent des supports de filtr^ 

pour mettre en oeuvre une methode de filtrage spa- 
10 tial par SELECTION-CGMBINAISON. 

les F1G.4A k 40 illustrent la combinaison des filtres 

des FIG.3A k 3D. 

les F1G.5A et 5D montrent des supports de filtres 
FMH. 

15 - les FtG.6A a 6C illustrent la combinaison des filtres 
des FIG.SA k 5D. 

les FIG. 7 A et 7B illustrent une methode de determi- 
nation de I'ecart type du bruit. 

20 1/ Dispositrf a ravons X 

En reference avec la FIG.1. un systeme de radio- 
graphie numerisee comprend une source 1 de rayons X 
; une table mobile 2 pour recevoir un patient ; un dispo- 

25 sitif intensificateur d'image 3 couple k un tube video 4, 
qui fournit des donnees d'images k un systeme 5 de 
traitement numerique d'image comprenant un micropro- 
cesseur Ce dernier comprend plusieurs sorties dont 
une sortie 6 couplee a un moniteur 7 pour la visualisa- 

30 tion des images radiographiques ou des images d'inten- 
site d'une zone observee du patient. 

Une image radiographique numerisee se presente 
sous forme d'une matrice bidimensionnelle de pixels. 
L'image radiographique numerisee peut contenir, par 

35 exemple. 512 x 512 ou 1024 x 1024 pixels codes sur 8 
bits ou 10 bits. Chaque pixel peut etre ainsi affecte d'un 
parmi 256 ou 1024 niveaux d'intensite. Par exemple les 
regions sombres ont un niveau de gris faible. et les 
regions claires de l'image ont un niveau de gris eleve. 

40 L'image num6ris6e peut etre obtenue en mode 
fluoroscopique de rayons X. L'invention peut trouver 
une application particulierement avantageuse dans le 
filtrage d'images angiographiques. Les prises de vue 
peuverrt §tre effectuees sous forme de sequence d'ima- 

45 ges a un taux de plusieurs images par second es. par 
exemple 30 images par seconde en mode fluoroscopi- 
que standard, ou 12,5 images par seconde en mode 
fluoroscopique pulse. 

L' elimination du bruit sans creation de TACHES DE 

50 BRUIT doit alors pouvoir etre faite sur chaque image de 
la sequence arrivant en temps reel. Le procede de fil- 
trage doit done §tre rapide, et ne pas detruire les 
details. 

Dans tous les cas. la presente invention ne tient 
55 pas compte de la methode par laquelle llmage numeri- 
que a ete obtenue, ni de la nature des objets qu'elle 
represente, mais concerne seulement le filtrage de 
cette image consideree seule et non en relation avec 
d'autres images de la sequence pour eliminer le bruit 
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sans generer de TACHES de bruit 
11/ Procede de fittraae spatial 

En reference avec la F1G.2, 'le prcx:ede de fiftrage 
spatial comprend la definition, dans une image J a trai- 
ler, d'un pixel appele pixel courant Po(3'o.yo) de coor- 
donnees Xo/q dans la nnatrice bidimensionnelle de 
rinnge J. Untensite lo(xo,yo) de ce pixel est supposee 
bruitee au depart. 

Pour un pixel courant Po(^.yo) d'intensite Iq 
(xoYo) dans llmage J, on def init. dans une fen§tre glis- 
sante WD(xo.yo) centree sur ce pixel courant, une serie 
formees d'un nombre N de firtres elementaires. 

Chacun des filtres de la serie est un filtre orients. 
Par filtre oriente. on entend qu'une dimension de son 
support est plus grande que I'autre et que cette grande 
dimension, ou longueur, est parailele ^ un axe appel6 
axe d'orientation, passant par le pixel courant 
Po{xo.yo). I'autre dimension, ou largeur etant perpendi- 
culaire ci cet axe. L'angle forme par I'axe d'orientation 
avec, par exemple. la direction OX des lignes de la 
matrice bidimensionnelle de I'lmage, est appele angle 
d'orientation O du tfltre eiementaire. 

Ainsi. la serie des fiHres elementaires lies au pixel 
courant dans ta fenStre WD(xo,yo) comprend plusieurs 
caract6ristiques parmi lesquelles. pour chaque filtre : 

son angle d'orientation 0, 
I'etendue de son support egale au nombre de pixels 
parallelemertt ^ I'axe d'orientation multiplie par le 
nombre de pixels perpendiculairement a cet axe, 

et pour la serie de filtres : 

le nombre N des filtres, 


Puis on realise une selection des pixels a nrrterieur des 
supports en fonction de leur intensite, pour selectionner 
" les pixels appartenarrt potentiellement au m§me objet 
que le pixel courant. Ensuite, on effectue une combinai- 

5 son lineaire ou moyenne des rntensites sur chacun des 
supports des filtres elementaires, et on combine les 
resultats de ces moyennes d'intensite, de maniere a 
foumir une intensite resuttante, qui est I'intensite fiitree 
k attribuer au pixel courant 

10 Dans ce procede de selection-combinaison, les 
supports des filtres incluent le pixel courant Po{xo.yo). 

En reference avec la FIG. 3 A, on a represerrte, ^ 
titre d'exemple, la d^inition d'une serie de N filtres 
orierTtes Fi, F2, F3, F4. dont, pour chaque filtre : 

15 

Tangle d'orientation ©t . ©2, ©3. ©4 est fixe, 
I'etendue du support est aleatoire, 
les supports de ces filtres sont sym^triques par rap- 
port au pixel courant Po(xo.yo). 

20 

et le nombre N de filtres est fixe et egai a 4. 
En reference avec la F1G.3B, on a represente. k 
titre d'exemple, la definition d'une serie de N filtres 
orientes F^. Fj. F3, F4,. dont, pour chaque filtre : 

25 

l'angle d'orientation Q^. ©2. ©3, ©4 est aleatoire, 
I'etendue du support est aleatoire. 
les supports sont dissymetriques par rapport au 
pixel courant Po{xo.yo). 'e nombre N de filtres 
30 est fixe et egal a 4. 

Le nombre N des filtres elementaires peut egale- 
ment etre aleatoire. 

En reference avec la FIG.3C, on a represente une 
35 serie de N = 3 filtres elementaires orientes FI, F2. F3, 
dont pour chaque filtre : 


Chaque filtre de la serie de filtres elementaires 
dans la fen^tre WD(xo,yo) liee au pixel courant 
Po(xo,yo). comprend au moins une caracteristique de 
valeur aleatoire. On peut ainsi definir, dans ia fendtre 
WD(xo.yo) centree sur le pixel courant Pofeyo). une 
serie de N filtres orientes. dont les angles d'orientation 
© sont aleatoires, les autres caracteristiques 6tant fixes. 
On peut aussi definir une serie de N filtres orientes, 
dont k la fois les angles d'orientation © sont aleatoires 
et les etendues des supports sont aleatoires. 

En outre, dans la serie, le nombre N de filtres peut 
Stre aleatoire. 

On decrit ci-apr^ plus predsement les caracteristi- 
ques et les proprietes du precede de fiitrage spatial 
dans lequel on definit differentes series de tels filtres 
elementaires aleatoires. 

11.1/ Procede de selection-combinaison 

Dans ce procede appele de selection-combinaison, 
on definit d'abord la structure des filtres elementaires 
orientes ayant une ou des caraaeristiques aleatoires. 


40 


l'angle d'orientation Q^. ©2. ©3 est aleatoire. 
I'etendue des supports est aleatoire, et les supports 
de ces filtres sont symetriques par rapport au pixel 
courant Po(xo.yo)' 


En reference avec la F1G.3C, on peut se fixer pour 
les filtres. un nombre ^^ax d'axes d'orientation poten- 

45 tieis fixes, regulierement repartis dans I'e^ce autour 
de Po(xo.yo)- choisissant un nombre N de filtre 
inferieur ou egal a N^^x. 'es filtres se repartissent de 
maniere aleatoire le long des N^iax d'orientation. 
Par exemple, on peut choisir N^^x = 8 ; dans ce cas. les 

50 8 axes d'orientation PqZ^. P0Z2. Pq^Z' - Pq^s poten- 
tiels possibles sont regulierement espacees de n/S ; et 
on peut choisir te nombre N de filtres fixe et egal k 3. 
Ces 3 filtres se repartissent alors de maniere aleatoire 
le long des 8 directions fixes possibles. 

55 Dans un autre exemple, on peut choisir le nombre 
N de filtres fixe et egal a 4. Ces filtres se repartissant de 
fagon aleatoire sur les 8 axes d'orientation. 

Dans ce type d'exemple iilustre par la FIG-3C, ou 
des axes d'orientation possibles sont fixes, le nombre N 
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desfiltres peut 6galement etre aleatoire. en restant infe- 
ri eur au nombre N^a, d'axes d'orientation. 

Dans chacune des F1G.3A et 38, I'etendue des 
supports des f iltres est de 1 x K pixels, ou K est le nom- 
tre de pixels dans le support-, ce nombre K etant al6a- 5 ' 
toire dans chacun de ces deux exemples. Les filtres 
representes a titre d'exemple ont done des supports 
monodimensionnels. Dans d'autres exemples non 
representes, les supports des filtres peuvent §tre bidi- 
mensionnels k condition qu'une des dimensions soit au io 
moins double de Pautre dimension, cette derni^re etant 
petite, par exemple 2 ou 3 pixels De preference, la 
fen§tre WD(xo.yo) est choisie de Tordre de 20 x 20 
pixels, ce qui fixe la grand e dimension des filtres k au 
plus 20 pixels dans ce cas. is 

Parmi les caracteristiques aleatoires potentielle- 
ment utilisables dans la definition de ia serie des N fil- 
tres elementaires orientes pour atteindre aux buts de 
llnvention. c'est-^-dire pour realiser un filtrage spatial 
du bruit sans generation de TACHES (PATCHES) et 20 
sans deterioration des details, contours et contraste, on 
propose dans I'ordre decroissant d'importance : 

1) Les support des N filtres a sont d'etendue 
aleatoire par exemple avec les orientations k 0n 25 
fixes et les autres caracteristiques fixes, comme 
illustre par la FiG.3A ; 

2) les supports des N filtres k sont d'etendue 
et d'orientation ©1 k ©2 aleatoires, les autres carac- 
teristiques etant fixes, comme illustre par la FIG. 38 30 

3) les supports des N filtres F^ a F^ sont d*etendue 
et d'orientation ©^ ^ ©^ aleatoire, et en outre le 
nombre N des filtres est aleatoire. comme illustre 
par la FiG.3C ; 35 

4) les supports son! en outre asymetriques, comme 
illustre par la FIG.3B ; 

5) un sous-echantillonnage dans la fenetre 
WD(xo.yo) peut dtre effectue prealablement autrai- 
tement de chaque pixel et le pixel de depart du 4o 
sous-echantilionnage peut etre aleatoire ; 

6) la moyenne MEAN pour obtenir soit les resultats 
individuels Ri. Rg, Rn des filtres. soit le resultat 
general R peut s'apptiquer sur des intensites pon- 
derees par des poids aleatoires ; k cet effet. on d^i- 45 
nit un coefficient aleatoire C lie au support, egal a 
rinverse de la somme des poids, compris errtre 0 et 

1 , affectes k chaque pixel dans le support ; chacun 
de ces poids peut etre egal k 1 : dans ce cas le 
coefficient est egal a Tinverse 1/K du nombre K de so 
pixels contenus dans le support ; le coefficient C de 
chaque filtre est pris en compte pour effectuer la 
combinaison lineaire ponderee des intensites. 

Par exemple le resultat Ri peut §tre obtenu en ss 
appliquant sur les differentes intensite 1^ (xq+u, yo+v) 
du filtre Fl des poids aleatoires et en faisant la 
moyenne MEAN de ces intensites ponderees. d'ou il 
resulte : 


R ^ = C x H [W ^(u.v) x I i(Xo+u. yQ+v)] 

Toutes les caracteristiques dtees plus haut peuvent 
etre rendues aleatoires. seules ou en combinaison, en 
sachant neanmoins que les caracteristiques les plus 
importantes sont. pour les supports de filtres I'etendue 
aleatoire et Torientation aleatoire et pour la combinaison 
le nombre de filtres N aleatoire. 

Exemple i : Le proced6 de selection-combinaison 
SIG-MEAN-MEAN 

Ce premier precede de SELECTION-COMBINAI- 
SON est d6signe ci-apr6s par le terme SIG-MEAN- 
MEAN pour simplifier Tecriture. 

Ce precede comprend d'abord, comme decrit ci- 
dessus, la definition des filtres orientes avec : 

leurs caracteristiques aleatoires, 
leurs caracteristiques fixes, 
leur nombre fixe ou aleatoire. 

Ce precede comprend ensuite la SELECTION dans 
chaque support de filtre d'un certain nombre de pixels 
en fonction de leur difference d'intensite vis-e-vis de 
I'intensite du pixel courant. 

Ce precede comprend ensuite la combinaison des 
filtres. 

En posant que le pixel courant Pq a pour coordon- 
nees (xo.yo) matrice bidimensionnelle de 

llmage J. et que le pixel P que Ton examine a pour coor- 
donnees (xq+u, yo+v), ou u et v sont des nombres 
entiers, la SELECTION est effectuee pour tout pixel P 
ayant une intensite l(xo+u. yo+v) telle que : 

IK^o+u. Yo+v) - io{Xo.yo)l < m X cTg 0) 

ou m est un facteur multiplicatif et ou oq est recart type 
du bruit dans t'image (noise standard deviation, en 
anglais). 

uecart type du bruit dans I'image peut etre estime 
par toute rnethode connue de retat de la technique, ou 
bien peut etre estime a priori. L'ecart type du bruit est 
mesure en niveaux de gris (ou intensite). II peut etre par 
exemple 5 ou 10 niveaux de gris. II est indispensable de 
realiser une estimation de l'ecart type du bruit ag pour 
mener a bien le precede selon I'invention. 

On rappelle ci-apres la definition de l'ecart type du 
bruit telle que connue de I'homme du metier 

uecart type du bruit note est l'ecart moyen que 
fait le bruit par rapport k sa moyenne. L'ecart type du 
bruit est la racine carree de la variance ctI du bruit 
dans I'image. L'ecart type du bruit peut etre obtenu en 
appliquant la rnethode fournie ci-apres a titre d'exemple. 

En reference avec le diagramme de blocs de la 
FIG. 7 A, cette rnethode comprend d'abord le traitement 
de I'image J par ligne. pixel par pixel, selon un balayage 
qui peut etre un balayage standard, depuis en haut k 
gauche de I'image jusqu'a en bas a droite de Hmage. 
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Dans une premiere etape symbolis^e par le bloc MEM 
100. les intensites l(x«y) des pixels sont memorisees en 
fonction de leurs coordonnees (x,y) dans la matrice bidi- 
mensionneile de I'innage. 

,Dans un^ seconde etape symbolisee par les blocs 5 
1 10A et 1 10B. et illustree par la FiG.78, en chaque pixel 
courarrt P(x,y)-de coordonnees x.y d'une ligne L de 
pixels dans Umage J. on calcule deux moyennes des 
intensites des pixels situ^ k gauche du pixel courant. si 
le sens de balayage de la ligne est vers la droite : io 

une premiere moyenne Ml des intensites sur un 
grand: support calcu! symbolise par le bloc 
1 10A;: c'est-a-dire une moyenne des intensites d'un 
grand nombre de pixels de la ligne situes a gauche is 
du pixel courarrt. ces pixels ayant deja ete balayes : 
par grand support, on entend environ 10 pixels, 
une seconde nioyenne M2 des intensites sur un 
petit support S2. calcul symbolise par le bloc 1 10B. 
c'est-a-dire une moyenne des. intensites d'un plus 20 
petit nombre-de pixels que cetui du grand support 
par exemple la mottie, soit 5 pixels, et a gauche du 
pixel courant 0(x,y) sur la ligne. 

Ces premiere et seconde moyennes d'intensites 25 
Ml. M2 peuvent §tre exactes ou ponderees. comme 
connu de I'homme du metier. Par ponderees, on entend 
qu'on affecte. k chaque intensite des pixels dont on fait 
la moyenne. un poids compris entre 0 et 1 . 

Dans une troisieme etape symbolisee par le bloc 30 
120A. la methode de determination de I'ecart type du 
bruit ag, comprend ensuite le calcul d'une valeur A1 
egale ^ la valeur absolue de la difference entre la pre- 
miere moyenne d'intensite M-i sur un grand support et la 
seconde moyenne d 'intensity M2 sur un petit support 35 
devaluation de cette difference : 

permet de rejeter les intervailes de pixels, sur ta tigne 40 
traitee qui contiennent des bords d'objet, c'est-^-dire 
dans lesqueis apparatt une brusque variation d'intensite 
due aux contours d*un objet : on declare qu'il y a un bord 
d'objet dans I'intervalle de pixels constitue par le grand 
support S-) , lorsque la seconde moyenne d'intensite M2 45 
sur le petit support est tres differente de la premiere 
moyenne d'intensite sur le grand support ; c'est-ci- 
dire quand la difference Ai calculee entre ces deux 
moyennes d'intensite excede une valeur estimee en 
niveaux de gris egale a environ 10 % des niveaux de so 
gris utilises pour construire Timage J numerique ; et, si 
un bord est detecte, on rejette I'intervalle constitue par 
les pixels du grand support Si . 

Dans une quatrieme etape symbolisee par le bloc 
120B, cette methode de determination de I'ecart type 55 
comprend en outre le calcul d'une valeur A2 egale a la 
valeur absolue de la difference entre la premiere 
moyenne d'intensite sur le grand support et I'lntensite 
du pixel courant l(x.y) 


A2 = |M, -l(x.y)| 

L' Evaluation de cette difference permet de rejeter 
les pixels aberrants. On dWare que le pixel courant a 
Une intensite abenante lorsque la difference d'intensite 
A2 llntensrte du pixel courant l(x,y) et la premiere 
moyenne d'intensite M^ sur un grand support 
depasse une valeur estimee en niveaux de gris comme 
d6finie plus haut. environ 10 % des niveaux de gris de 
Hmage J. 

Dans une etape ulterieure symtDOlis^e par le bloc 
de TEST 130. cette methode comprend I'examen des 
deux conditions precedentes : 

si le pixel courant n'a pas une intensite aberrante, 
si Untervalle de pixels correspondant au grand sup- 
port ne contient pas de bord d'objet. 

Dans le cas ou ces deux conditions sont conjoint e- 
ment realisees. alors on declare que A2 qui est la valeur 
absolue de la difference d'intensite entre le pixel cou- 
rant l(x.y) et la premiere moyenne d'intensite sur le 
grand support S^, est une mesure elemental re a(x,y) 
retenue pour revaluation de I'ecart type du bruit 

La suite des operations decrites ci-dessus est repe- 
tee pour tous les pixels k traiter de I'lmage J. Puis une 
moyenne general e M de toutes les mesures elemental - 
res est effectuee dans le bloc 140 et est retenue comme 
valeur de I'ecart type du bruit ag. 

Dans la relation (1). le facteur m multiplicatif de 
I'ecart type du bruit ag peut etre egal a 2 ou 3. Dans le 
cas ou le bruit est gaussien et ou un filtre elementaire 
est applique sur une zone de I'inoage ou i'intensite est 
sensiblemerrt uniforme, si m = 2, alors environ 3/4 des 
pixels sont sElectionnes, et si m = 3 alors pratiquement 
tous les pixels sont selectionnes. 

On examine alors. en appliquant la relation (1 ), tous 
les pixels situes k I'interieur de chaque support du filtre 
elementaire oriente. 

Loperation de SELECTION decrite ci-dessus, qui 
prend en compte I'ecart type du bruit (lire SIGMA B) 
pour selectionner des pixels d'intensite proche de celle 
du pixel courant est designee ci-apr6s par le terme BIG 
pour simplifier I'ecriture. 

En reference avec la FIG. 3D. on a montre dans une 
image J. une feneti-e WD(xo.yo) centree sur le pixel 
courant Po(Xo.yo). des supports de filtres Fi. F2. F3, F4, 
monodimensionnels orientes dans diverses directions 
du plan de i'image. autour du pixel courant Po(xotyo)' 
on a mis en relief les pixels selectionnes par SIG en les 
entourant d'un cercle. 

A Tissue de I'operation de SELECTION SIG, une 
combinaison lineaire ou moyenne des intensites des 
pixels selectionnes est effectuee. Cette operation de 
moyennage est designee ci-apres par le terme MEAN. 

En reference avec la FiG.4A qui est un diagramme 
du filtrage par SELECTION-COMBINAISON SIG- 
MEAN-MEAN d'une serie de filtres telle qu'illustree par 
la FIG.3D, le filtre F, effectue la selection des pixels 
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symbotisees par le bloc SIG 10.1. et calcule dans le 
bloc MEAN 11.1 une premiere moyenne des diver- 
ses intensites 

li' (xq-kj. Vo+v) des pixels selectionnes dans F, avec : 
'r, =MEAN(i/(Xo-KJ.yo+v)] 

tel que 

II/Cxq-hJ. yo+v)-lo(=^o*yo)l <^^<^b {^) 

(xq+u, Vo+v) decrivant le support oriente du fiitre F^. 

Puis, le fiitre Fg calcule dans le bloc MEAN 11.2 
une seconde moyenne R2 des diverses intensites 
l2'{xo+u. yo+v) des pixels selectionnes dans un bloc 
SIG 10.2 avec : 

R2=MEAN [la'Cxo+u. yo+v)] 

ou I'intensite des pixels selectionnes obeit k la relation 
(1). (xq+u. yo+v) decrivant !e support oriente du fiitre 

On opere de meme jusqu'au fiitre F^ compris (SIG 
10.N. MEAN 11. N) ; 

R^ = MEAN [lN'(Xo+u.yo+v)l 

ou ies diverses intensites 1n'(xo+^. Yo+v) des pixels 
selectionnes obeissent a la relation (1). (xq+u, yo+v) 
decrivant le support orients du fiitre F^. 

Lorsque toutes les moyennes Hi ^ R^ ont et6 obte- 
nues par les filtres orientes F^ k F^. une moyenne 
MEAN est effectuee dans MEAN 12 sur lesdttes moyen- 
nes Ri a Rn. ce qui fournit I'intensite filtree RoC^^.yo) 
resuitante a appliquer finatement au pixel courant 
Po(xo.yo) 

Ro(^o.yo) = MEAN (R1.R2-RN] 

Exemple II : Le precede de selection comblnaison 
SIG-MEAN-MED 

Ce second proced6 de selection-coinbinaison est 
designe ci-apr§s pour la simplidte d'^criture par le 
terme SIG-MEAN-MED. A Tissue de I'operation de defi- 
nition de la serie de N filtres ^lementaires orientes avec 


leurs caracteristiques aleatoires, 
ieurs caracteristiques fixes, 
leur nombre fixe ou aleatoire, rcper'ation de selec- 
tion des pixels dans chaque support de fiitre est 
effectuee comme dans I'exemple I par le crit^e (1) 
portant sur Tintensite des pixels et prenant en 
compte Tecart type du bruit cg. 

Ensuite. en reference avec la F1G.4B, Tobtention 
des resultats R^ a R^ particuliers des filtres F, a F^ est 
effectuee par un moyennage MEAN sur les diverses 
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intensites respectives I/Cxq+u. yo+v) a In'(xo+o, yo+v) 
des filtres k F^. puis les intensites R^ k Rn resultant 
de ces moyennages sont combinees par un fiitre 
median designe par MED 13, qui foumit I'intensite resul- 
5 . tante Ro{xo.yo) telle que : 

Ro(Xo.yo) = MED [R1.R2... Rn) 

Un fiitre median est un fiitre qui effectue le classe- 
10 ment des intensites qui lui sont presentees, avec un 
numero d'ordre, de la plus petite a la plus grande inten- 
site, et qui fournit comme resultat. Ilntensite dont la 
valeur est celle du numero situe au milieu de I'ensemble 
des numeros, ou si les numeros ne sont pas impairs, la 
15 moyenne dintensite coaespondant aux deux numeros 
du centre. 

L'avantage de la methode SIG-MEAN-MED est que 
le fiitre median utilise a des proprietes d'injection du 
bruit haute frequence qui le rendent specialement inte- 

20 ressant pour eviter i'effet de TACHES (PATCHES) dans 
rimage. En effet le fiitre median effectue la selection 
d'un des resultats, entre R1 et RN qui est la valeur 
situee au milieu de ces valeurs d'intensite, au lieu d'une 
moyenne de ces intensites. Le resultat Ro{xo.yo) Qui 

25 est obtenu avec le MEDIAN est une intensite plus forte- 
ment decorrelee. d'un pixel courant traite, au pixel cou- 
rant suivant traites, que I'intensite resultant de la 
moyenne MEAN 12 utilisee dans I'exemple I decrit pre- 
cedemment. 

30 

Exemple III : Le proced^ de selection comblnaison 
KNN 

Ce troisieme precede des selection combinaison 
35 est designe ci-apr6s, pour la simplicite d'ecrrture par te 
terme KNN. 

Dans cet exemple, une autre methode de SELEC- 
TION en fonction de I'intensite est appliquee au lieu du 
crit^re (1) qui tient compte de I'ecart type du bruit 03*. 
40 Cette autre methode de selection consiste a retenir tous 
les pixels dont I'intensite differe le moins de celle du 
pixel courant. A cet effet, on calcule les differences 
d'intensite 

45 A(Xo.yo) = IKxq+u, yo+v) - 1 o(Xo.yo)l (2) 

entre les pixels dans chacun des filtres et I'intensite du 
pixel courant. 

Dans cette methode appelee des K PLUS PRO- 
50 CHES VOISINS et designee par le terme KNN, on 
selectionne en fait un nombre fixe K de pixels ayant 
I'intensite la plus proche de celle du pixel courant 

Toutes les autres caracteristiques et etapes du pro- 
cede KNN peuvent etre identiques a celles de I'exemple 
55 I ou I'exemple II, et sont illustrees k cet effet par les 
FIG.4C et 4D. 

Dans ces FIGS, les blocs KNN 20.1 a KNN 20.N 
symbolisent les etapes de selection des pixels, et les 
blocs MEAN 21.1^ MEAN 21 .N symbolisent la realisa- 
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tion des moyennes d'intensite des pixels selectionnes 
sur les fiitres. 

Les blocs MEAN 22.1. et MED 22.2. effectuent res- 
pectivenrient une moyenne ou un filtrage median des 
resyltats particuliers des fOtres elementaires pour four- s 
ntr le resultat final Ro(^yo)- 

11.2/ Le precede FMH 

Dans les etapes du precede decrit ci-apres on d^i- w 
nit une serie de N fiitres lineaires elementaires moyen- 
neurs. k reponse impulsionnelle. dont les reporrses sort 
en outre combinees par des fiitres medians, ce qui 
revient k designer cette classe de fiitres par FMH (en 
anglais : FIR HYBRID MEDIAN FILTERS), Cette combi- is 
naison a une structure d'arbre qui comprend les fiitres 
lineaires mcyenneurs (FIR) et les fiitres medians 
(MED). 

En reference avec les FIG.5A k 50, on definit une 
fenetre WD(xo,yo) bidimensionnelle, centree sur le 20 
pixel courant Po(xo.yo) ^® coordonnees (xQ.yo) dans la 
matrice bidimensionnelle de I'lmage J a traiter. Dans 
cette fengtre WD(xo,yo). on definit des supports de fii- 
tres selon les lignes des pixels, et les colonnes de 
pixels, et les deux diagonal es entre les lignes et tes 2S 
colonnes. Le pixel central PoCxo'Vo) ^® fendtre 
WD(xo,yo), qui est le pixel courant. est considere k 
part. 

En reference avec ta FIG. 5 A, on a represente a titre 
d'exemple une serie de N fiitres orientes Hi. H2. V^, V2. 30 
D^ 1 . D-12, D21 . D22 dont pour chaque f iltre : 


En reference avec la FIG.5B, une serie de fiitres du 
meme type est representee avec les caracteristiques 
suivantes, par rapport a la FIG.5A : 


les supports des fiitres qui sont disposes selon la 
meme direction sont dissymetriques de mani^re 
aleatoire. 

' En reference avec les F1G.5C et 50, au lieu de con- 
aderer un nonnbre de fiitres disposes selon tous les 
angles potentieis d'orientation e^. 02. 03. ©4. seule- 
ment deux axes sont consider^es : les axes horizontal 
et vertical, sur la FiG.SC ; les axes formes par les direc- 
tions des deux diagonales sur la F1G.5D. D*ou, sur la 
FIG.5C. la serie de fiitres comprend seulement les sup- 
ports Hi. H2. Vi. V2 alignes sur ces axes horizontal et 
vertical ; et sur la FiG,5D, la serie de fiitres comprend 
seulement les supports D^, D^g. D21 et D22 alignes sur 
ces axes en diagonale. Sur ces F1G.5C et 5D k titre 
d'exemple. les supports ont une etendue aleatoire et 
sont symetriques par rapport au pixel courant 

En reference avec les FIG.6A a 60, les fOtres ayant 
les supports de caracteristiques aleatoires decrrts en 
relation avec les F1G.5A ^ 5D effectuent d'abord la com- 
binaison lin^aire ou moyenne MEAN des intensites des 
pixels sur leur support. 

On a appel6 IV1. 1'v2. ^'hi- 1'h2 '^s differentes inten- 
sites des pixels situes sur les supports Vi. V2. H^ et H2 
respectivenrient. 

On a appele I'dh, t'Di2. l'o2i» '*d22 differentes 
intensites des pixels situes sur les supports D11, D12, 
D21. D22 respectivement. 

Les r^ultats des moyennes effectuees. dans les 
blocs MEAN31.1 a MEAN 31.8. par les fiitres H^ H2, 
Vi. V2, D11, D-j2. D21, D22 sont des intensites notees 
'hi- 'h2» 'vi. 'v2' 'oil' 'di2' 'o2i 'd22 respectivement. 

En reference avec la FIG.6A, la combinaison en 
arbre des resultats des fiitres moyenneurs MEAN pre- 
cedents comprend 2 niveaux induant I'introduction 
dans un filtre median MED 32,5 a 5 valeurs, des inten- 
sites moyennes lyi. iv2. 'hi* 'h2 sur les axes horizontal 
et vertical plus TirTtensite de Iq (Xq, yo) clu pixel central ; 
et en parallel e, Tintroduction dans un filtre median MED 
32.6 a 5 valeurs, des intensites moyennes Iqh. Idi2' 
'd21 . b22 sur les deux directions diagonales plus I'inten- 
site \q{xq, yo) du pixel central. Ces deux fiitres medians 
fburnissent respectivement les intensites resultantes 
R1 et R2. terminant le premier niveau de I'arbre. 

Les deux intensites R1 et R2 calculees au premier 
niveau de I'arbre sont introduites, avec Tintensite lo{xo. 
Yo) du pixel central dans un nouveau filtre median MED 
33. a 3 valeurs, qui fournit I'intensite resultante 
RQ(xo,yo) ^ appliquer comme valeur d'intensite fiitree 
sur te pixel courant Pq (xq. yo)- 

En reference avec la FIG.6B. la combinaison com- 
prend egalement 2 niveaux, mais eile n'est appliquee 
que sur les moyennes d'intensite lvi. Iv2- 'hi. 1h2 relati- 
ves aux seules directions verticales et horizontales res- 
pectivement. Cette combinaison comprend 
introduction des moyennes d'intensite lyi, Iv2 P'*^^ 
I'intensite lo(^.yo) bruitee du pixel central dans un filtre 
median MED 32.1 a 3 valeurs pour fournir une intensite 
resultante R1 : et I'introduction des moyennes d'inten- 


I'orientation est fixe consideree par example par 
rapport k I'axe OX, selon la verticale 0^. Thorizon- 
tale ©2 (angle "O" represente) et les diago'nales 3S 
©3, ©4 passant par le pixel courant centre du carre 
formant la fen§tre, 

r etendue du support est aleatoire : ces supports 
sont monodimensionnels. c'est-^-dire ont une lar- 
geur de 1 pixel et une longueur qui est limitee par le 40 
cadre de la fenetre WD (xo.yo). 
les supports des fiitres Hi, H2 disposes selon la 
meme direction, par exenrple horizontale. sent 
symetriques. II en est de mdme des supports de fii- 
tres et V2 disposes selon la direction verticale. 45 
des supports de fiitres D-i-, et D-]2 disposes selon la 
premiere diagonale. des supports de fiitres D^2> 
D22 disposes selon la secorKle diagonale. 
le nombre N des fiitres est fixe ; ici ce nombre est 8. 
et ces supports de fiitres sont disposes selon tous so 
les angles possibles d'orientation 61 a ©4 multiples 
de Tz/4 dans I'espace, 

dans la serie de fiitres. on incorpore le pixel central 
Po{xo,yo) de la fenetre et son intensite est conside- 
ree a part 55 
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site Ihi. Ih2 P^^ilele dans un fdtre median MED 32.2 
a 3 valeurs pour foumir une intensite resultante R2. 

Les Intenstes R1 et R2 caJcuiees au preniier 
niveau de I'arbre sont introduites, avec I'lntensite lo(xo. 
yo) du pix^! central dans un nouveau filtre median MED 5 
33.1, k 3 valeurs, qui fournit rintensit6 resultante 
RQ(xo,yo) ^ appliquer comme valeur dlntensite filtree 
sur le pixel courant Pq {xq. yo)- 

On peut dans une combinaison non representee^ 
remplacer les valeurs de moyennes d'intensite Ivi, Iv2- 
'hi< Ih2 ^^"s FIG.66 par les moyennes obtertues 
dans les memes conditions sur les deux diagonales, 
soit respectivement Iqh. 'di2> '021. 'd22' 9^^^ former 
une autre combinaison de f litres k deux niveaux sur 
seulement deux axes qui sont celles des deux diagona- is 
les. 

En reference avec la FIG.6C. la combinaison en 
arbre comprend 3 niveaux. Les deux premiers niveaux 
sont constitues comme decrit en relation avec la FIG.6B 
de la combinaison des filtres sur les deux axes horizon- 20 
tal et vertical et donne par MED 33.1 un resultat R3 : et 
en parall^le comme decrrt plus haut. de la combinaison 
des filtres sur les deux axes diagonaux et donne par 
MED 33.2 un resultat R4. 

Le premier resultat R3 et le second resultat R4 sont 25 
ensuite combines avec I'lntensite bruit^e Iq (xq. yo) 
pixel courarrt, dans un filtre median MED 34 ^ 3 valeurs. 
pour donner en sortie de ce troisidme niveau, I'intensite 
r§suitante Ro(^.yo) ^ appliquer comme valeur d'inten- 
site filtree sur le pixel courant Pq (xq. y©)- 

Lintensite d'un pixel courant est filtree par la serie 
de filtres'elementaires ai§atoires prealablement definis, 
selon les etapes de selection et de combinaison propo- 
sees, puis les Intensites des autres pixels desires de 
llmage sont filtrees, jusqu'ci ce que tous les pixels desi- 35 
res de I'image aient subi le traitement sefon le precede 
propose, c'est-a-dire. en operant le filtrage de I'intensite 
de chaque pixel desire par une serie de filtres elemen- 
taires dont les caracteristiques qui sont fixes ou qui sont 
al6atoires ont ^te prealablement definies. ainsi que le ^ 
mode de realisation des combinaisons d'intensite. 

Revendications 

1 . Proced6 de filtrage spatial du bruit dans une image 45 
(J) ayant la forme d'une matrice bidimensionnelle 
de pixels [P(x.y)] ayant des valeurs d'intensite 
[l(x,y)] numerisees bnjitees, ce precede compre- 
nant des etapes de : 

50 

a) selection, dans une fen&tre glissante 
[WD(xo, yo)] centree sur un pixel courant 
(Po(xo. yo)] d'u" nombre (N ^ 0) de filtres ele- 
mentaires (Fi. F2...- F(sj ; V^, \/2, H-i. H2, On, 
* D12. D21, D22] orientes aleatoires ayant des 55 
supports avec un axe d'orientation [PqZ^, 

P0Z2 Pq^] passant par le pixel courant et 

des caracteristiques dont une ou plusieurs ont 
des valeurs aleatoires. 
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b) combinaison des intensites des pixels sur 
lesdits supports individuellement. foumissant 
une intensite el6mentaire, et combinaison des 
intensites elementaires, fournissant une inten- 
. sit6 filtree (Ro(Xo. yj] ^ appfiquer au pixel cou- 
rant. 

2. Precede selon la re^endication 1, selon lequel les 
supports des filtres ont potentiellement. par mi les 
caracteristiques dont une ou plusieurs ont des 
valeurs aleatoires : 

leur nombre (N) dans la serie, 

un angle d'orientation . ©g.... ©jsj entre I'axe 

d'orientation [Po^i . ^0^2 Pq^nJ et un axe 

de reference (OX) defini dans la matrice bidi- 
rectionnelle, 

une etendue definie en nombre de pixels dont 
la plus grande dimension est paralieie a I'axe 
d'orientation, 

une asymetrie par rapport au pixel courant, 
un coefficient Ii6 a des coefficients de pondera- 
tion appliques aux intensites des pixels. 

3. Precede selon Tune des revendications 1 ou 2. 
selon lequel un sous-echantillonnage des pixels de 
ia fenetre est potentiellement pratique avec un pixel 
de depart de ce sous-echantillonnage localise de 
maniere aleatoire. 

4. Precede selon I'une des revendications 1 k 3, seion 
lequel lesdits supports de filtre sont monodimen- 
siennels avec les pixels atignes sur les axes 
d'orientation. ou sont bidimensionnels avec la 
dimension du support parallel e k I'axe d'orientation 
au moins double de ia dimension transverse. 

5. Precede selon I'une des revendications 1 a 4, selon 
lequel : 

les angles d'orientation potentiels [©i, ©2.-.. 
©(sj] sont fixes et reguiierement repartis dans 
I'espace, 

le nombre (N). fixe ou aleatoire, des filtres ele- 
mentaires aleatoires est inferieur ou egal au 
nombre (N^nax) angles d'orientation. 
et les axes d'orientation des supports de ces fil- 
tres elementaires sont disposes selon ces 
angles d'orientation de maniere fixe ou alea- 
toire. 

6. Precede selon la revendication 5, selon lequel ; 

les supports des filtres elementaires aleatoi- 
res incluent le pixel courant 
ce precede comprenant en outre, entre retape a) et 
retape b). une etape de : 

a') selection (SIG. KNN) des pixels en fonction 
de leur intensite sur chaque support de filtre 


9 


EP 0 755 023 A1 ^ 18 


elementaire. selon un entire qui retient les 
pixeis dintensite proche de celle du pixel cou- 
rant pour estimer leur appartenance k un 
meme objet de Timaga 

5 

7. Prcxede selon la revendication 6, selon lequel. pour 
la selection des pixels, te crrtere retient les pixels 
qui ont une intensite [I(xo-hj. yo+v)] telle que la 
valeur absdue de la difference entre ladtte intensite 

et llntensite (lo(xo. yo)l clu pixel courant est tnfe- io 
rieure k un multiple (m) de l*ecart type du bruit (ctb) . 
dans I'inriage (J). 

8. Precede selon ia revendication 6, selon lequel : 

pour ia selection des pixels, le critere retient is 
les pixels en un nombre predetermine (K) dans le 
support- 

9. Precede selon I'une des revendications 7 ou 8, 
selon lequel la combinaison des intensiles des 20 
pixels sur les supports est Iin6aire et la combinai- 
son des intensites elementaires est effectuee par 

un filtre moyenneur (MEAN). 

10. Precede selon Tune des revendications 7 ou 8, 25 
selon lequel (a combinaison des intensites des 
pixels sur les supports est lineaire et la combinai- 
son des intensites elenientaires est effectuee par 

un filtre median (MED). 

30 

11. Precede selon ia revendication 5. selon lequel : 

les axes d'orientatiorl fixes sent au nombre de 
quatre, disposes selon la ligne, la coionne, et 
les deux diagonales passant par le pixel cou- 35 
rant, dans la fen§tre, 

les supports sent positionnes pour exclure le 
pixel courant. etant disposes de part et d'autre 
du pixel courant. sur la ligne (H^ H2), sur la 
coionne (Vi. V2). sur la premiere diagonaie 40 

D12) et sur la seconde diagonaie {D21. 
D22) passant par le pixel courant, 
le pixel courant, central dans la fenetre [Wq 
(xq. yo)]. est consider^ comme un support par- 
ticulier, 45 

et la combinaison des intensites elemental- 
res prend en compte Tintensite du pixel courant. 

12. Precede selon la revendication 11, selon lequel ia so 
combinaison des intensites des pixels sur les sup- 
ports est lineaire. et la combinaison des intensites 
elementaires est effectuee par des filtres medians 
sur deux ou trois niveaux, prenant en compte 
I'intensite du pixel courant. ss 

13. Precede selon Tune des revendications 1 a 12. 
selon lequel les etapes definies pour un pixel k trai- 
ter sent repetees pour tous les pixels a traiter dans 


llmage. 

14. Disposrtif pour Taide ^ Tobservation d'une image 
medical e. connprenant : 

un syst^me pour foumir des donnees d'une 
image sous forme de matrice bidimensionnelle 
de pixels ayant des valeurs d'intensrte numeri- 
sees, 

un syst^me d'affichage pour I'affichage de 
rimage. 

un systeme de traitement dimages ayant 
acc^s aux donnees d'image. et au systeme 
d'affichage, 

ce systeme de traitement d'image etant arrange 
pour effectuer des etapes selon le precede con- 
forme a Tune des revendications 1^13. 
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